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We consider multitype branching processes arising in the study of random laminations of 
the disk. We classify these processes according to their subcritical or supercritical behavior and 
provide Kolmogorov-type estimates in the critical case corresponding to the random recursive 
lamination process of [1]. The proofs use the infinite dimensional Perron- Frobenius theory and 
quasi-stationary distributions. 
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1 Introduction 

In this note we are interested in multitype branching processes that arise in the study of random 
recursive laminations. In order to introduce and motivate our results, let us briefly recall the 
basic construction of [I]. Consider a sequence U\, V\, U2, V2, ■ ■ ■ of independent random variables, 
which are uniformly distributed over the unit circle Si . We then construct inductively a sequence 
Li, L2, ■ ■ ■ of random closed subsets of the closed unit disk D. To start with, L\ is set to be the 
(Euclidean) chord [£/iVi] with endpoints U\ and V\. Then at step n + 1, we consider two cases. 
Either the chord [£7 n +iV„+i] intersects L n , and we put L n +i = L n . Or the chord [J7 n +lKi+l] 
does not intersect L n , and we put L n+ \ = L n U [U n+ iV n+ i]. Thus, for every integer n ^ 1, L n 
is a disjoint union of random chords. See Fig. 1. 







Figure 1: An illustration of the process creating the sequence (L n ) n ^i. We use hyperbolic chords 
rather than Euclidean chords for aesthetic reasons. 

A fragment of L n is a connected component of D\L n . These fragments have a natural 
genealogy that we now describe. The first fragment, D, is represented by 0. Then the first 
chord [LqVi] splits D into two fragments, which are viewed as the offspring of 0. We then 
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order these fragments in a random way: With probability 1/2, the first child of 0, which is 
represented by 0, corresponds to the largest fragment and the second child, which is represented 
by 1, corresponds to the other fragment. With probability 1/2 we do the contrary. We then 
iterate this device (see Fig. 2) so that each fragment appearing during the splitting process is 
labeled by an element of the infinite binary tree 

T 2 =\J{0,l} n , where {O,l}° = {0}. 

n^O 

If F is a fragment, we call end of F, any connected component of F fl Si. For convenience, 
the full disk D is viewed as a fragment with end. Consequently, we can associate to any 
u G T2 a label £(u) that corresponds to the number of ends of the corresponding fragment in 
the above process. Lemma 5.5 of [1] then entails that this random labeling of T2 is described by 
the following branching mechanism: For any u £ T2 labeled m ^ 0, choose m\ G {0, 1, . . . , m} 
uniformly at random and assign the values 1 + m\ and 1 + m — m\ to the two children of u. 
This is the multitype branching process we will be interested in. See Fig. 2. 




Figure 2: On the left-hand side, the first 7 chords of the splitting process. On the right-hand side, the 
associated branching process corresponding to the number of ends of the fragments at their creations. 
Notice that we split the fragments according to the order of appearance of the chords, thus the binary 
tree on the right-hand side seems stretched. 

We can also define a random labeling by using the above branching mechanism but starting 
with a value a ^ at the root of T2, the probability distribution of this process will be denoted 
P a and its relative expectation E a . A ray is an infinite geodesic path u = (m, U2, . . .) € {0, 1} N 
starting from the root in T2. For any ray u = (u\, . . . , u n , . . .) or any word of finite length 
u = (ui, . . . , u n ), we denote by [u]j or [u]i the word (ui, . . . , U{) for 1 ^ i ^ n, and [u]q = 0. 

Theorem ([ , Lemma 5.5]). Almost surely, there exists no ray u along which all the labels 
starting from 4 are bigger than or equal to A, 

P 4 (3u G {0,1} N : £([u]i) > 4,Vt > 0) = 0. 

The starting label 4 does not play any special role and can be replaced by any value bigger 
than 4. This theorem was proved and used in [1] to study certain properties of the random closed 
subset Lqo = UL n , and in particular to prove that it is almost surely a maximal lamination 
(roughly speaking that the complement of is made of disjoint triangles), see [1 , Proposition 



2 



5.4] . One of the purposes of this note is to provide quantitative estimates related to this theorem. 
Specifically let 

G n = {ue {0, l} n : £([u]i) > 4, Vi G {0, 1, . . . , n}} 

be the set of paths in T2 joining the root to the level n along which the labels are bigger than 
or equal to 4. 

Theorem 1.1. The expected number of paths starting from the root and reaching level n along 
which the labels starting from 4 are bigger than or equal to 4 satisfies 

E 4 [#G n ] -> -^i-. (1) 

Furthermore, there exist two constants < ci < C2 < 00 suc/i i/ia£ i/ie probability that G n 7^ 
satisfies 

01 < P 4 (G n /0) < C2 . (2) 

Remark 1.2. These estimates are reminiscent of the critical case for Galton-Watson processes 
with finite variance a 2 < 00. Indeed if H n denotes the number of vertices at height n in such a 
process then E [H n ] = 1 and Kolmogorov's estimate [ ] implies that P (H n 7^ 0) ~ -§-. 

The proof of Theorem 1.1 relies on identifying the quasi-stationary distribution of the labels 
along a fixed ray conditioned to stay bigger than or equal to 4. This is done in Section 2. 
In Section 3, we also study analogues of this branching random walk on the fc-ary tree, for 
k ^ 3, coming from a natural generalization of the process (L n ) n ^o where we replace chords by 
triangles, squares... see Fig. 3. 






Figure 3: Extension of the process (L n ) n ^i where we throw triangles or squares instead of chords. 

We prove in these cases that there is no critical value playing the role of 4 in the binary case. 

Acknowledgments. The first author thanks Microsoft Research and the University of 
Washington, where most of this work was done, for their hospitality. We are also grateful to 
Jean-Frangois Le Gall for precious comments and suggestions on a first version of this note. 



2 The critical case 
2.1 A martingale 

Fix an arbitrary ray uo in T2, for example uo = (0, 0, 0, 0, 0, . . .) and define X n = £([uo] n ) for 
n ^ 0, so that X n is the value at the n-th vertex on the fixed ray uq of the T2-indexed walk I 
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starting from xq 4 at the root. Then (X n ) n ^o is a homogeneous Markov chain with transition 
probabilities given by 

P2(x,y) = — — rli^+i. 

We first recall some results derived in [I]. If J- n is the canonical filtration of (X n ) n ^Q then 
a straightforward calculation leads to E XQ [X n +i \ F n \ = 1 + X n /2, hence the process M n = 
2 n {X n — 2) is a martingale starting from xo — 2. For i ^ 1, we let Tj be the stopping time 
Ti = inf{n ^ : X n = i}, and T = T\ A T2 A T3. By the stopping theorem applied to the 
martingale (M n ) n ^Q, we obtain for every n ^ 0, 

x - 2 = ~E X0 [M nAT ] = E X0 [-2 Tl l{ Tl=T ^ n }] + +E a;o [2 T3 l{ T3=T ^ n }] + ~E X0 [2 n (X n - 2)l{ T>n }]. 

One can easily check from the transition kernel of the Markov chain (X n ) n ^Q that for every 
i ^ 1, P Xo [Ti = T = i] = P XQ [T2 = T = i] = P Xo [T^ = T = i] . Hence, the equality in the last 
display becomes 

x - 2 = F, XQ [2 n (X n - 2)l| T>n }], 

or equivalently 

x -2 = 2 n P X0 [T>n] B X0 [X n -2 \ T > n}. (3) 

Our strategy here is to compute the stationary distribution of X n conditionally on the non 
extinction event {T > n}, in order to prove the convergence of ~E4[X n \ T > n] and finally to 
get asymptotics for P^T > n}. Before any calculation, we make a couple of simple remarks. 
Obviously E X() [X n - 2 | T > n] ^ 2, and thus we get 2 n P Xo (T > n) sC S^. Since there are 
exactly 2 n paths joining the root of T2 to the level n, we deduce that the number #G n of 
paths joining to the level n along which the labels are bigger than or equal to 4 satisfies 

E, [#G n ] < (4) 

Notice that a simple argument shows that if 4 ^ xo ^ xi then the chain X n starting from xo 
and the chain X' n starting from x\ can be coupled in such a way that X n ^ X' n for all n ^ 0. 



2.2 The quasi- stationary distribution 

We consider the substochastic matrix of the Markov chain X n killed when it reaches 1,2 or 3: 
This is the matrix (P2(x, y)) x ^A given by 

h(x,y) = ^-j-lj^+i- 

We will show that P2 is a 2-recurrent positive matrix, in the sense of [3, Lemma 1]. For 
that purpose we seek left and right non-negative eigenvectors of P2 for the eigenvalue 1/2. In 
other words we look for two sequences {g(x)) x ^4, and (f(x)) x ^4 of non-negative real numbers 
such that /(4) = g(4) = 1 (normalization) and for every x ^ 4 

a{x) = 2 s ^g(y)P 2 {y,x) 

f{x) = 2^P 2 (x,y)f(y) 

y>4 



2 £ 



j / =(x-l)V4 
n x+1 



a(y) 
y + i ! 



X + 



i £/(*)• 



(5) 
(6) 



y=4 
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We start with the left eigenvector g. From (5), we get g(5) = <?(4) = 0, and g(i) — g{% + 1) 
\g(i — 1) for % 5. Letting 



z i+i 



the last observations lead to the following differential equation for G 

2G(z) = z- 1 (z-l)G'{z) + z 3 , 

with the condition G(z) = z 5 /5 + o(z 5 ). A simple computation yields G{z) = 3/4exp(2z)(z — 
l) 2 + (z 3 /2 + 3z 2 /4 — 3/4). After normalization, the generating function Gi/ 2 (z) = J2i>4 9l/2('i')z i 
of the unique probability distribution g-yM which is a left eigenvector for the eigenvalue 1/2 is 
given by 

Gi/ 2 (*) = |(exp(2z)(z-l) + z + l) ! 

that is 

This left eigenvector is called the quasi- stationary distribution of X n conditioned on non- 
extinction. For the right eigenvector /, a similar approach using generating functions is possible, 
but it is also easy to check by induction that 

/(*) = -jj— lt>4, 

satisfies (6). Hence the condition (iii) of Lemma 1 in [3] is fulfilled and the substochastic 
matrix P 2 is 2-recurrent positive. For every x ^ 4, set q n {x) = ~Pi(X n = x \ T > n) = P^(T > 
n)~ 1 P 2 n v (x) where v stands for the "vector" (/Uj)i^4 with = 1 and t> j = if i ^ 5. Theorem 
3.1 of [3] then implies that 

q n {x) — > Si/20)- (7) 



Unfortunately this convergence does not immediately imply that E4 [X n \ T > n] — > E [X] 
where X is distributed according to <?x/2- But this will follow from the next proposition. 

Proposition 2.1. For every n the sequence I — - I is decreasing. 



51/2O) 



x>4 



Proof. By induction on n ^ 0. For n = the statement is true. Suppose it holds for n 0. By 
the definition of q n +i, for x 4 we have 

<7„+i(x) = P 4 (X n+ i = x I T > n + 1) 

= p 4(r + = *n+i = s, T>n) 

P 4 (T > n) v 
P 4 (T >n + l)^z + l lJ 

We need to verify that, for every x ^ 4, we have q n +l( x )9i/2( x + 1) ^ q n+ \(x + l)g\/2{ x ) or 
equivalently, using (8) and (5) with g = gi/ 2 , that 

( y gl/2(*A / yv gnQg) I > ( ffl/2(f) J / yv gnfcA 

\2;^2;V 4 Z + 1 / \z>(a:-l)V4 * + / W(a;-l)V4 2 + 1 I WiV4 Z + V 



For x = 4 this inequality holds. Otherwise, if x > 4, we have to prove that 



, . #1/20) ^ g n (ar) 

7 „(x-i) 2^ ~7qrf > 51/2(^-1) 2^ j^y- ( 9 ) 



Set ^ = ^^—^ to simplify notation. The induction hypothesis guarantees that q n (z) ^ 
A x gi/ 2 (z) for every z ^ x, and therefore 

gn(g) < A ffi/20) 

^ z + 1 ^ * ^ z+1 ' 



This gives the bound (9) and completes the proof of the proposition. □ 

By Proposition 2.1 we have for every x ^ 1, ^li^ <; c, where C = sup n>0 < 

oo by (7). This allows us to apply dominated convergence to get 



e 2 + 3 

E 4 [X n |T > n] = ^xq n {x) -—^^xg l/2 {x) = G' 1/2 (l) = ■ 



n— >oo 

a;^4 a;>4 



Using (3) we then conclude that 



2 n P 4 [T>n] > - 9 -. (10) 



. - e 2 . 

2.3 Proof of Theorem 1.1 

(n) 

We first introduce some notation. We denote the tree T2 truncated at level n by Tj . For 
every u = (ui, . . . , u n ) G {0, 1}™, and every j G {0, 1, ... , n}, recall that [u]j = (m, . . . , Uj), 
and if j ^ 1, also set [u]* = (ui, . . . , Uj-i, 1 — Uj). We say that j G {0, 1, . . . ,n — 1} is a 
te/f tern (resp. turn) of ii if = (resp. uj + i = 1). A down step of u is a time 
j G {0, 1, . . . , n — 1} such that 

e([u} j )-e([ u ] j+l )^2. 

Note that if j is a down step of u then ^([it]* +1 ) = 2 + £([u]j) — £([u]j+i) ^ 4. The set of all 
j G {0, 1, . . . , n — 1} that are left turns, resp. right turns, resp. down steps, of u is denoted by 
L(u), resp. R(it), respD(u). We endow T2 with the lexicographical order ^ , and say that a path 
u G {0, l} n is on the left (resp. right) of v G {0, l} n if u ^ v (resp. v ^ u). A vertex of {0, l} n 

In) 

will be identified with the path it defines in T 2 ■ If u, v G T2 we let aAti be the last common 
ancestor of u and u. 

Proof of Theorem 1.1. Lower bound. We use a second moment method. Recall that 
G n = {u£ {0, l} n : *([ti]i) > 4, Vi G {0, 1, . . . , n}} 

(n) 

is the set of all paths in T 2 from the root to the level n along which the labels are bigger than 
or equal to 4. A path in G n is called "good". Using (10), we can compute the expected number 
of good paths and get 

E 4 [#G n ] = 2™P 4 [T>n] — ► -^i-, 

n— >oa e — 1 

as n — y 00, which proves the convergence (1) in the theorem. For u G G n and j G {0, 1, . . . , re}, 
we let Right (u,j) be the set of all good paths to the right of u that diverge from u at level j, 

Right (it, j) = {v G G n : u <v and u A v = [u]j}. 
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In particular, if j is a right turn for u, that is itj+i = 1, then Right (it, j) = 0. Furthermore 
Right (it, n) = {u}. Let us fix a path it G {0, l} n , and condition on u G G n and on the labels 
along u. Let j G {0, 1,2,..., n}. Note that the first vertex of a path in Right (it, j) that is not 
an ancestor of u is and its label is 2 + — £([u]j + \), so if we want Right (u, j) to be 

non-empty, the time j must be a down step of u. If j is a left turn and a down step for u, the 



subtree {w G Tg : it; A [u]* = [it]?} on the right of [it] 3 - is a copy of T^" J ±; , whose labeling 
starts at ^([«]* +1 ). Hence thanks to (4) we get 



(n-j-l) 



E 4 [# Rights, j) | « G G n , (^([u]i)) <i<»] < 



[u]j)-£([u] j+1 ) 



Since the labels along the ancestral line of u cannot increase by more that one at each step, if 
u G G n we have Y^7=o I ^(N«+i) ~~ ^(Mi) I IjeD(u) ^ n - Combining these inequalities, we obtain 



E 4 



^#Right(u,j) ueG n , (£{[u]i)) ^ 

3=0 



n 



We can now bound E 4 [#G^] from above: 



E 4 [#G^] < 2E 4 



ue{o,i} n u<v 



= 2 P4(nGG„)E 4 

«e{o,i} n 

= 2 P4(nGG„)E 4 

■ue{o,i} n 
^ n. 



veG„ 



u G G n 



^#Right(n,j) 

3=0 



u G G n 



fill 



The lower bound of Theorem 1.1 directly follows from the second moment method : Using (1) 
and (11) we get the existence of c\ > such that 



P(#G w >0)^ _ E[#G "L 



E[(#G r 



n 



(12) 



Upper Bound. We will first provide estimates on the number of down steps of a fixed path 
u G {0, l} n . Recall that L(u), R(it) and D(u) respectively denote the left turns, right turns, and 
down steps times of u. 

Lemma 2.2. There exists a constant C3 > such that, for every n ^ and every uq G {0, l} n 

P 4 (u G G n , #D(u ) < c 3 n) < C 3 1 2~ n exp(-c 3 n). 

Proof. We use the notation of Section 2.1. For any set A C {0,1,..., n — 1} and m G 
{0, 1, . . . ,n - #A}, with the notation JV^ = #{j G {0, 1, ... ,n - 1}\A : = 5} we have 
from [ , formula (27)] 



P[X j+1 > (Xj-l) V4, Vj G {0, 1, . . . , n — l}\vl , JV^ = m] < 



1 \ m / 3 \ n—m—#A 



We will first obtain crude estimates for iV^. Note that N A ^ iV^f and that sup^ 4 P2(i, 5) = |, 
so that for any 5 C {0, 1, . . . , n} we have 



P (Xj = 5 , Vi G B) < 5" 
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By summing this bound over all choices of B with #B ^ m we get P (N® ^ m) ^ 2 n 5 m for 
every m E {0, 1, . . . n}. Let ki E (0, 1/2) and «2 £ (0, 1) such that K\ + k 2 < 1. We have 

P(u E G n , #D(n ) ^ Kin) 
^ P (n E G n , # D(uo) < Km , JV® < K 2 n) + P (iVjf > K 2 n) 

P [X j+1 > (Xj — 1) V 4, Vj E {0, 1, . . . , n — 1}\A ; < K 2 n] + P « ^ rc 2 n) 



Ac{0,l,..,n-1} 



^ (L^nJ+1) 



< n 



n 

\_KlTl 



Ac{0,l,...,n-1} 
[K 2 n\ 



n— L^inJ 



+ P (N% ^ K 2 n) 



n— \_K\n\ 



(13) 



Notice that for every A > 1 we can choose /«i > small enough so that nL^ n <j yl n for 
large enough. Furthermore 



n— [nin\ 



2 n 2 



g\ n-LKinJ-L«2nJ 

7 J 



and by choosing k\ even smaller if necessary we can ensure that the right hand side of (13) is 
bounded by Cg 2~ n exp(— c^n) for some C3 > 0. □ 



We use the last lemma to deduce that 

nP 4 (3tt £G„, # D(n) < c 3 n) < 



?i 



C3 



exp(-c 3 n) 



0. 



(14) 



We now argue on the event El = {3u E G n , #(D(u) n L(u)) c^n/2}. On this event there 
exists a path u E G n with at least c^n/2 down steps which are also left turns. Conditionally on 
this event we consider the left-most path P of G n satisfying these properties, that is 

P = mm{u£ G n , #(D(u) n L(tt)) > c 3 n/2}. 
A moment's thought shows that conditionally on P and on the values of the labels along the 

(n) 

ancestral line of P, the subtrees of T 2 hanging on the right-hand side of P, that are the 
offsprings of the points [P]j+i for j E L(P), are independent and distributed as labeled trees 
started at £([P]* j+1 ). 

Hence conditionally on P and on the labels ((£([P]j),0 ^ i ^ n), for any j E L(P) n D(P) 
the expected number of paths belonging to the set Right (P, j) (defined in the proof of the lower 
bound) is 



E, 



#Right(P,j) P, (£([P]i)W<n 



= 2^- 1 P, ([P] * +i) (T>n-i-l) 
^ 2 n -J- 1 P 4 (T>n-j-l) 

> K 3 > 0, 



(15) 



where K3 is a positive constant independent of n whose existence follows from (10). Thus we 
have 



E4[#G n I El] 



E, 



E 



^#Right(P,j) P, (£([P],)) 
_i=o 

> K 3 E 4 [#(D(P)nL(P))|E L ]. 



. C3K3 

n. 



(16) 



Since P±(El) ^ E4[#G n ]/E4[#G n | El] we can use (1) to obtain ~P±(El) ^ n^/n for some 
constant K4 > 0. By a symmetry argument, the same bound holds for the event Er = {3u G 
G n , #(D(u) n R(u)) ^ C3?i/2}. Since {G n / 0} is the union of the events Er, El and 
{3u G G n , # D(-u) ^ C371}, we easily deduce the upper bound of the theorem from the previous 
considerations and (14). □ 



3 Extensions 

Fix k ^ 2. We can extend the recursive construction presented in the introduction by throwing 
polygons instead of chords: This will yield an analogue of the multitype branching process on the 
full fc-ary tree. Formally if X\ . . . . . Xfc £1X6 k (distinct) points of Si we denote by Pol(a;i, . . ■ ,x k ) 
the convex closure of {x±, . . . ,x k } in D. Let (Uij : 1 ^ j ^ k , i ^ 1) be independent random 
variables that are uniformly distributed over Si. We construct inductively a sequence L\,L k , ■ ■ ■ 
of random closed subsets of the closed unit disk D. To start with, L\ is Pol^i^, . . . , U\^). 
Then at step n+ 1, we consider two cases. Either the polygon P„+i := Pol(C/" n +i t i , . . . , f7 n +i,jfc) 
intersects L*, and we put = L^. Or the polygon P n +i does not intersect L*, and we 

put L^ + i = U Pfc. Thus, for every integer n ^ 1, L\ is a disjoint union of random /c-gons. 
In a way very similar to what we did in the introduction we can identify the genealogy of the 
fragments appearing during this process with the complete fc-ary tree 

T fc = U{0,1,...,A;-1}*, where {O,l,...,£:-l} o = {0}. 

Then the number of ends of the fragments created during this process gives a labeling l k of T fc 
whose distribution can be described inductively by the following branching mechanism (this is 
an easy extension of [1, Lemma 5.5]): For u G Tfc with label m^Owe choose a decomposition 
m = mi + ^2 + . . . + rrik with mi, ni2, • • • , G {0, 1, . . . , m}, uniformly at random among all 
("^-l" 1 ) P ossl ble choices, and we assign the labels mi + l,rrt2 + 1, . . . , + 1 to the children 
of 0. Again the distribution of the labeling i k of obtained if we use the above branching 
mechanism but started from a ^ at the root will be denoted by P a and its expectation by E Q . 
We use the same notation as in the binary case and are interested in a similar question: For 
which a ^ does there exist with positive probability a ray u such that £ fc ([u]j) ^ a for every 
i 0? Specifically, the value a is called subcritical for the process (£ k (u),u GTj.) when there 
exists a constant c > such that 

P a (3u G {0,1,. . ., k - l} n : l k {[u]i) > a, Vi G {0, 1, . . . ,n}) < exp(-cn). 

It is called supercritical when there exists a constant c > such that we have both 

f P a (3uG{0,l,...,fc-l} N :^([u] 4 ) ^a, Vi G {0, 1, . . .}) > c, 

\ E fl [#{ue{0,l,... I fe-l} n :^(Hi)>a, Vi G {0, 1, . . . , n}}] ^ exp(cn). 

Note that a deterministic argument shows that if k 2 and a = 2, there always exists a ray 
with labels greater than or equal to 2, also when k = 2 and a = 3 there exists a ray with labels 
greater than 3. The case k = 2 and a = 4 has been treated in our main theorem. We have the 
following classification of all remaining cases: 

Theorem 3.1. We have the following properties for the process & , 

• for k = 2 and a ^ 5 the process is subcritical, 

• for k = 3 the process is subcritical for a ^ 4, and supercritical for a = 3, 

• for k ^ 4 and a ^ 3 i/ie process is subcritical. 
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Proof. Supercritical Case k = 3 and a = 3. We will prove that for k = 3 and a = 3, the 
process is supercritical. Similarly as in Section 2.1 we consider the tree- indexed process £ 3 on a 
fixed ray of T3, say {0, 0, 0, . . .}. Then the process Y n given by the n-th value of £^ started from 
3 along this ray is a homogeneous Markov chain with transition matrix given by 



2(x + 2- y) 
(x + l)(x + 2) 



We introduce the stopping times Tj = inf{n ^ 0, Y n = i} for i = 1, 2 and set T = Ti A T2. We 
consider a modification of the process that we denote Y" n , which has the same transition 
probabilities as Y n on {1,2,3,4} , but the transition between 4 and 5 for Y n is replaced by a 
transition from 4 to 4 for Y n . Thus we have Y n ^ 4 and an easy coupling argument shows that 
we can construct Y n and Y n simultaneously in such a way that Y n ^ Y n for all n ^ 0. Hence 
we have the following stochastic inequality 

T < T, 

with an obvious notation for T. To evaluate T we consider the subprocess Y nA f which is again 
a Markov chain whose transition matrix restricted to {3,4} is 

1/5 1/10 
1/5 1/5 

The largest eigenvalue A max of this matrix is greater than 0.34, which implies that 

P 3 T > n ^ P 3 (T > n) > K 5 (0.34) n , 

for some constant K5 > independent of n. It follows that the expected number of paths starting 
at the root of T3 that have labels greater than or equal to 3 up to level n, which is 3 n P (T > n), 
eventually becomes strictly greater that 1: There exists no ^ 1 such that Pa(T > n) > 3~ n for 
n ^ riQ. A simple coupling argument shows that the process £ 3 started from a 3 stochastically 
dominates the process £ 3 started from 3. Consequently, if we restrict our attention to the levels 
that are multiple of no and declare that v is a descendant of u if along the geodesic between u and 
v the labels of £ s are larger than 3, then this restriction stochastically dominates a supercritical 
Galton- Watson process. Hence the value 3 is supercritical for £ s . 
Subcritical Case k = 3 and a = 4. As in the binary case we let 



2(x + 2-y) 
(x + l)(x + 2) 



be the substochastic matrix of the process Y n started at 4 and killed when it hits 1,2 or 3. We 
will construct a positive vector (h(x)) x -^4 such that J2 X < 00 anc ^ 

h-P 3 < Xh, (17) 

for some positive A < 1/3, where we use the notation h ■ P${y) = ^2 x h{x)P'i(x,y). This will 
imply that 

P(T>n) < ^|^A", 

where T is the first hitting time of {1,2,3} by the process Y n started at 4. The subcriticality 
of the case k = 3 and a = 4 follows from the preceding bound since there are 3 n paths up to 
level n and A < 1/3. To show the existence of a positive vector x satisfying (17) we begin by 
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studying the largest eigenvalue of a finite approximation of the infinite matrix P3. To be precise 
let Pg 30 ^ = (Ps(i, i))4<ij/<30- A numerical computation with Maple® gives 



At, 



:= max { Eigenvalues (P 3 (30) )} ~ 0.248376642883065 < 1/3. 



The vector (h(x)) x ^4 is then constructed as follows. Let (h(x))^ x <^3o be an eigenvector asso- 
ciated with the largest eigenvalue A max of P^ 3 °\ such that min^a^o h(x) = /i(30) = 1. Note 
that the vector h can be chosen to have positive coordinates by the Perron-Frobenius theorem 
and it is easy to verify that x —> h{x) is decreasing. For x ^ 31 we then let 



h(x) = 13 



x-30 



30! x 2 



xi 



We now verify that this vector satisfies (17) with A slightly greater than A m ax* 

Suppose first 

that y G {4, . . . , 30}. In this case ^4<a;<30 ^( x )^3( x > v) equals A max /i(y) by definition, whereas 
the contribution of X^3i h(x)Ps(x, v) is l ess than Yl x >3i h(x) < 0.014, thus 

h-P 3 (y) ^ 0.263%). (18) 

Now, if y ^ 31 we have 

2 /or>l\ 2 



£ h(x)P 3 (x,y) *S 13^ 30 lu-'Psiy -l,j)(r^) + P 3 (y, y) + £ 13^' 



< 10,-30^0^7 2 ^ 2 4 v ,/ y\ V\ 

\yl) \l3y(y + l) + (» + l)(y + 2) + ^ VCy + OV ^ 

< 0.3-%) 

which proves (17). 

Other critical cases. The other critical cases are treated in the same way. We only provide 
the reader with the numerical values of the maximal eigenvalues of the truncated substochastic 
matrices that are very good approximations of the maximal eigenvalues of the infinite matrices, 

max{eigenvalues(P 2 (i,j))5<i,i<3Q} ^ 0.433040861268365 < 1/2, 

max{eigenvalues(Pi(«,i))3<i,i<3o} 0.231280689028977 < 1/4. 

□ 
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